£tA L b-4-~j /■? 



h £7k’'4-i7a2; 

H tn Vfv ^ £1 ’T.J 




5Sir 



r- 141 

Copy 

RM L54I19 


?-'V. 


NACA 


RESEARCH MEMORANDUM 


FLUTTER INVESTIGATION OF A TRUE -SPEED DYNAMIC MODEL WITH 
VARIOUS TIP-TANK |:ONFIGURATIONS 

h 

Sewall, Robert W.^Herr, and William B. Igoe 

Langley Aeronautical Laboratory 
Langley Field, Va. 

'OS ara ®- 

CLASSIFIED DOCUMENT 

This material contains information affecting the I^tional Defense of the United States within the meaning 
: the espionage laws, Title 18, U.S.C., Secs. 793 and 794, the transmission or revelation of which in any 
iamner to an unauthorized person is prohibited by law. 

TIONAL ADVISORY COMMITTEE 



FOR AERONAUTICS 

WASHINGTON 

March 30, 1955 


NACA RM L54I19 


CONTIDENTIAL 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

FLUTTER INVESTIGATION OF A TRUE-SPEED DYNAMIC MODEL WITH 
VARIOUS TIP-TANK CONFIGURATIONS 
By John L. Sewall, Robert W. Herr, and William B. Igoe 

SUMMARY 


A l/6-scale dynamic model, equipped with wing tip tanks and repre- 
sentative of current \mswept fighter-type airplanes of conventional plan 
form, was tested for flutter. The model was dynamically scaled to flut- 
ter at the same speed as a typical full-scale conf igviration . The param- 
eters important to flutter were satisfactorily approximated by the use 
of a spar -balsa segment-type wing construction. 

A flutter-arresting device was located within each tip tank which 
would provide for a sudden shift in the tip-tank center of gravity. This 
device was very effective in stopping symmetric flutter for cases where 
the tip-tank center of gravity was shifted forward of the elastic axis 
of the wing. 

The flutter investigation was concerned with the effects of external 
stores whose inertial and geometric properties were systematically varied. 
In general, the results obtained concur with those found in previous 
researches in that symmetric flutter speeds were increased for tip-tank 
centers of gravity forward of the wing elastic axis, and tended to decrease 
as the ratio of uncoupled wing bending frequency to torsion frequency 
approached vinity. Increasing the voltune of the tip tanks by a factor of 5 
tended to decrease the flutter speed. The addition of a horizontal fin 
at the rear end of each tank may have been beneficial for the smaller tank 
and sli g htly detrimental for the larger tank. 

Flutter speeds were calculated by means of a conventional Rayleigh- 
Ritz type of flutter analysis employing uncoupled modes . The effects of 
various assumptions in the calculations involving structural damping, 
higher structural modes. Compressible aerodynamic coefficients, and sting- 
mount flexibility were explored. The calculations were, in general, exces- 
sively conservative for tip-tank centers of gravity near the elastic axis; 
however, the introduction of structural damping materially improved the 
agreement between experimental and calciilated flutter speeds for ratios 
of imcoupled wing bending frequency to torsion frequency near unity. The 
introduction of sting-movint flexibility showed that this could have a 
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strong effect on flutter speed particularly in those cases where the 
natural frequency of the sting system was near the flutter frequency, 
and f\irther work in this regard is warranted in future flutter studies 
of sting -mounted models. Also indicated as desirable is the need for 
continued research on the oscillating aerodynamic coefficients for flutter 
calciilations on wings equipped with large external stores or nacelles . 


INTRODUCTION 


Within the past several years a number of cases of flutter are known 
to have occurred on fighter and light bomber airplanes (at least in their 
development stages) carrying large external stores. Considerable interest 
has, therefore, developed in flutter problems for wings carrying such 
stores . References 1 to 10 are representative of both the theoretical 
and experimental investigations of problems of this kind. The calcula- 
tions of flutter instabilities for these cases have, in general, been 
unsatisfactory where two-dimensional aerodynamic forces and moments 
together with a limited number of structural modes have been employed, 
and where aerodynamic forces and moments on the tip tank have been ignored. 
Because of weaknesses inherent in many of the more complete analyses avail- 
able at present, greater emphasis has come to be placed on wind-tunnel 
testing of dynamically scaled models with external stores . 

Although valuable experience has been obtained from flutter tests 
with scaled speed models (see, for example, refs. 11 to 15) > these tests 
cannot provide needed information on the effects of Mach number . The 
extrapolations from low-speed model results to full-scale high Mach num- 
ber conditions have proved troublesome and ambiguous, particularly where 
the margins of safety from flutter may have been small. In order to 
circumvent this difficulty, it is necessary to scale the model so that 
its Mach niimber at flutter will be the same as that of the prototype, or, 
in other words, to design and test a true-speed model. Accordingly, two 
principal objectives of the investigation reported in this paper were 
(l) to determine the feasibility of designing and building, a true-speed 
dynamically similar model with tip tanks to meet a given set of elastic 
and inertial specifications, and (2) to locate the experimental flutter 
boundaries of this model for a variety of tip-tank conditions. Another 
objective was to evaluate a flutter-arresting device involving a rapid 
shift of the tip-tank center of gravity and having application to full- 
scale-flight flutter-testing techniques . ' 

For the design of the model, in addition to keeping its flutter 
speed the same as that of the prototype, the principles followed in the 
scaling required that nondimensional parameters important in flutter, 
such as mass ratio, frequency ratio, and reduced frequency, remain the 
same on the model as on the full-scale counterpart. No attempt was made 
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to scale the stiffnesses in the tip-tank or wing -tip-tank attachment . 
Recent investigations reported in references and 15 have dealt with 
the Inertial and damping properties of sloshing fuel; however, in the 
present investigation fuel was simulated by means of solid weights, and 
no effort was made to account for the effects of fuel sloshing. 

The mnin effort Of the wind-tunnel test program was directed toward 
location of the flutter boundaries for a wide range of tip-tank inertial 
conditions corresponding to ratios of uncoupled wing bending frequency 
to torsion frequency near unity. Variations in these tip-tank conditions 
were accomplished by the use of fixed, solid weights. Also investigated 
were the aerodynamic effects of tip-tank size and the effects of a hori- 
zontal fin mounted at the rear portion of each tip tank. In order to 
study asymmetric flutter, the model was provided with one degree of free- 
dom in roll. 

The results of the flutter exi>eriments are compared with the resilLts 
of flutter calculations performed by utilizing a Rayleigh-Ritz type of 
analysis employing two-dimensional aerodynamic forces and moments together 
with a combination of uncoupled vibration modes to represent the flutter 
mode . No aerodynamic forces and moments were assumed to act on the tip 
tank. 


The scaled wing stiffnesses and inertial properties specified in 
the design of the model were chosen to be fairly representative of present- 
day fighter-type airplanes . Scaled wing stiffnesses and inertial prop- 
erties of an actual fighter airplane (the Lockheed F-9^B) are included 
as a comparison. Also included for purposes of comparison are the results 
of a flight flutter experience on a related fighter airplane (the Lockheed 
F- 8 OA) involving special tip-tank inertial loading conditions . 


SIMBOIS 


a nondimens ional wing elastic axis location relative to half- 

chord, positive for elastic axis rearward of midchord 

b half -chord of wing, ft 

bp flange width on wing spar, in. (see fig. 9) 

bj, reference half -chord of wing (taken at the station at three- 

fourths of the semispan) , ft 

web spacing on wing spar, in. (see fig. 9) 


CONFIDENTIAL 



CONFIDENTIAL 


NACA RM L54I19 


root chord of wing along fvaselage center line, in. (see fig. 5(a)) 

tip chord of wing at y = = 36, in. (see fig. 3(a)) 

distance between wing elastic axis and center of gravity, in.; 
positive for center of gravity rearward 

distance between wing elastic axis and tip-tank center of 
gravity, in.; positive for center of gravity rearward 

frequency, oi/2it, cps 

natural first bending frequency (experimental) , cps 
natural second bending frequency (experimental), cps 
flutter frequency, cps 

uncoupled first bending frequency ( calculated) , cps 

uncoupled second bending frequency ( calculated) , cps 

uncoupled frequency of effective sting-fuselage combination, 
cps (see appendix) 

natural first torsional frequency (experimental) , cps 

uncoupled first torsional frequency (calculated), cps 

uncoupled second torsional frequency (calculated) , cps 

coefficient of structural damping based on logarithmic decrement 
/ 1 ^ 

g = — logp — where a_ is the initial anrolitude and an 
rw ® an ° . 

the nth amplitude in n cycles of freely decaying oscillation] 
flange spacing on wing spar, in. (see fig. 9) 
reduced frequency, bto/v 


CONFIDENTIAL 



NACA RM L54I19 


CONFIDENTIAL 


5 . 


kg spring constant of equivalent sting support system, Ib/in. 

(see the appendix) 

I general dimension of length (see text on selection of scale 

factors) 

wing semispan from center line of fuselage to center of wing- 
tip-tajik attachments , in . 

m mass of wing per unit length, slugs/ft 

rg distance from effective pitch axis of rotation of sting to wing 

elastic axis at fuselage center line (see the appendix) 

ro, nondimens ional radius of gyration of wing section about elastic 

axis, |/lo,/mb2 

tp flange thickness on wing spar, in. (see fig. 9) 

web thickness on wing spar, in. (see fig. 9) 

v^ flutter speed, fps 

w weight of wing per unit length, Ib/in. 

Xq, nondilnens ional center -of -gravity location of wing section 

relative to half -chord, e/b; positive for center of gravity 
rearward of wing elastic axis 

y spanwise coordinate of wing as measured from fuselage center 

line , in . 

El wing bending stiffness, lb-in 

F + iG oscillating aerodynamic functions for two-dimensional incom- 

pressible flow (see ref. l6) 

GJ wing torsional stiffness, lb-in. ^ 


If mass rolling moment of inertia of one-half fuselage about 

^ fuselage center line, ft-lb-sec^ 

If mass pitching moment of inertia of one-half fuselage and effec- 

® tiye sting combination about effective axis of rotation in 

sting, ft-lb-sec^ (see the appendix) 
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Itr 

la 

lot 

M 

R 

W 

Wt 


W, 


w 


•asym 


■nsym 


mass rolling moment of inertia of tip tank about wing tip, 
ft-lb-sec^ 

mass pitching moment of inertia of wing per unit length about 
wing elastic axis, lb-sec^ 

mass pitching moment of inertia of tip tank about wing elastic 
axis, in-lb-sec^ 

Mach number 

Reynolds number 

weight, lb . 

weight of tip tank, lb 

weight of each wing panel, lb 

nondimens ional span coordinate used in antisymmetric flutter 
calculations , y/ Iw 

nondimens ional span coordinate used in symmetric flutter 

y - 3 


calculations , 


^w “ 3 


wing mass ratio, jtpb^/m 
scale factor for length, or 

tip-tank weight ratio, 

density of test medium, slugs/cu ft 
angular frequency, 2nf, radians/sec 
Subscripts : 

f flutter except as noted 

h bending degree of freedom 

r reference station at three-fourths semispan 

t tip tank 
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w wing 

F full-scale airplane except as noted 

M model 

S sting support system 

W spar web 

a torsional degree of freedom 

asym antisymmetric fuselage boundary condition (fioselage free to 

roll) considered in flutter calculations 

sym symmetric fuselage bovindary condition (cantilever and fuselage 

pitch) considered in flutter calculations 

DESCRIPTION OF MODEL 

The model (see fig. 1) was sting -mounted in the Langley l6-foot 
transonic tunnel (fig. 1(a)) and was representative of current unswept- 
wing fighter-type airplanes equipped with wing -tip fuel tanks . The aspect 
ratio of the wings without tip tanks was 6.0, the taper ratio was 
and the airfoil shape was an NACA 65AOI3 section. The wing was untwisted, 
and the 52-percent -chord line was normal to the longitudinal center line 

1 o 

of the airplane; although the elastic axis was swept back 3^ , the wing 
was considered vmswept . 

The model was tested with two different -size tip tanks (shown in 
fig. 2) which were mounted so that the longitudinal axis of the tank fell 
in the plane of the wing. Tests were conducted with and without hori- 
zontal fins located at the rear of the tip tank as shown in figure 2. 

Some wing details may be seen in figure 3* Each wing panel was mounted 
as a cantilever (see fig. 3(a)) to a cylindrical body simulating the 
fuselage as a rigid structiire that provided freedom of the model in roll. 


Selection of Scale Factors 

Scaling of the model was based on the simplified model construction 
discussed in reference 17. The qiiantities scaled were linear dimensions, 
mass, mass unbalance, moments of inertia, frequencies, and stiffnesses 
of the prototype . The wind tunnel ixsed had an octagonal, slotted throat 
and was capable of attaining Mach nximbers up to I.05. 
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The following table lists the flutter parameters considered signif- 
icant in this investigation in terms of the geometric scale factor A. 
The model was chosen to be l/6 the scale of a representative fighter- 
type airplane, and, therefore. 


X = = 1 or X = = 1 

bp 6 Zp D 

where the subscripts M and F refer to model and fxill scale, respec- 
tively. The choice of this particular value of A was influenced by 
the desire to obtain Reynolds numbers that would be as realistic as 
possible with respect to the full-scale airplane within the restriction 
imposed by available wind-tvinnel size . The density factor Pj^y^Pp was 

assumed to be unity. 


Parameter 

Symbol 

Scale factor 

Mass per unit length .... 

mp 

It 

Mass 

°M^m 

a 5 = -1- 


mpZp 

216 

Mass moment of inertia 


1. 1 

per unit length 

A^ - 



1296 

Mass moment of inertia . . . 

^“M^M 




7776 

Frequency 

“M 

1= 6 


CJCp 

A 

Bending stiffness 

(El)i^ 

A^ = — i— 


(ei)f 

1296 

Torsional stiffness .... 

(gj)m 

it 

H 


(GJ)ji 

1296 


(Sa)^ 


Structural damping 

( ^h) p ( So,) p 

1 

Mach number 

% 

1 


Mp 
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It should be noted that = Mp is not precisely true because of 

differences -which may exist between the stagnation temperature in the 
tunnel and that in free air for an airplane. No attempt was made to 
design a given value of structural damping (gjj^,ga) into the model. 

Accurate measurements of the bending and torsional damping coefficients, 
as determined from logarithmic decrements during the test program, were 
difficult. The coefficients measured varied from 0.01 to 0.10. The mean 
value appeared to be approximately 0.035 both bending and torsion 
degrees of freedom. This compares favorably with a structural damping 
value of 0 .044 obtained from a ctirrent fighter airplane, 'the 
Lockheed F-94B. 


Model-Construction Details 

Wings.- The solid c\xrves sho-wn in figures 4 to 8 give the scaled 
stiffness and inertial properties typical of current fighter-type air- 
planes, and the model was designed to meet these specifications. 

The wings were of spar-segment construction as shown in figure 5- 
The spar, the main structural member, was essentially a built-up box 
section made from thin sheet steel and welded at all four corners as 
sho-vm in section A-A of figure 3(a) . This type of spar construction was 
chosen because of the relative ease with which the bending and torsional 
stiffnesses could be varied essentially independent of one another. 
Dimensions of the spar are given in figure 9. 

The -wing plan form and airfoil shape were obtained by assembling a 
series of adjacent wing segments to the spar in the manner shown in fig- 
ure 3(b) for a typical segment. Ballast weights of the kind indicated 
in the figure were added to each segment to adjust the wing mass, mass 
mbalance, and mass moment of inertia to the design values. The segments 
were built-up hollow sections made from sheet balsa and were coated inside 
and out with fiber glass and Paraplex to provide extra rigidity for the 
segment to withstand better the forces encountered during violent flutter. 
Each segment was fastened to the spar flanges by means of four screws, 
each screw bearing on a small metal -backed strip of hardwood glued to 
the inside of the segment. 

The measiored inertial properties of the completed model are Indicated 
in figures 4, 5, and 6, and the spanwise distributions of wing stiffnesses 
in bending and torsion are given in figures 7 8, respectively. The 

attachment of the balsa segments to the spar caused no measurable Increases 
in bending and torsional stiffnesses . The gaps between the balsa segments 
were covered with stiffened paper for most of the flutter runs, and later 
this material was replaced -with silk for the remaining runs. Although 
the stiffness contribution of the silk was more noticeable than that of 
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the paper^ the measured stiffnesses of the assembled model were not 
significantly altered because of the lose of either material. However, 
the measvired stiffnesses which are presented in figures 7 and 8 were 
obtained prior to the flutter tests without either paper or silk seals 
attached . 

Stiffness measurements made at different times during the construc- 
tion of the model and the flutter test program revealed some weakening 
of the wing due to flutter violence; however, on the basis of frequency 
checks made during the flutter program, it appeared that such changes 
had little effect on the vibration characteristics of the wing. The 
dashed-line cvirves in figures 7 and 8 were considered most representative 
of the stiffness properties of the wing and were vised in the flutter 
calcvilations . 

Scaled stiffness and inertia parameters of a current fighter air- 
plane, the Lockheed F-9^B, are included in figures 4 to 8 for comparison 
purposes. For the inertial parameters (figs. 4 to 6) , these curves cor- 
respond to fvill fuel load in the integral wing fuel tanks . 

Tip tanks .- Flgvire 2 shows the two tip tanks that were used in this 
investigation. The smaller of these two tanks (hereinafter referred to 
as tank A) was a geometrically scaled reproduction of a full-sized Fletcher 
230-gallon tip tank, whereas the larger tank (tank B) was designed to 
have the same shape but three times the Internal volume of the smaller 
tip tank. No attempt was made to scale the stiffness of the tip tank or 
its attachment to the wing. 

Construction of the model tip tanks consisted of a duralumin cylin- 
drical center section with hollow plastic nose and tail cones made of 
fiber glass and Paraplex. Tank B was assembled by replacing the smaller 
nose and tail cones with larger ones and mounting a thin duralumin cover 
(shown removed in fig. 2) over the center section. For each tank, two 
tail cones were provided, one with a horizontal fin and one without a 
fin, in order to explore the effects of the fin on flutter . For tank A 
the fin was scaled geometrically from the full-sized configuration. Each 
tip tank housed a flutter -arresting device which is described in detail 
later in the paper. This device involved the use of one of three differ- 
ent shifting weights to achieve a rapid change in the mass imbalance of 
the tip tank. Table I gives the inertial properties of the tank for each 
of these weight conditions with no additional weights added. 

The geometrical properties of both tip tanks and fins are also given 
in table I. The fairing at the juncture of the wing and tank A c^ be 
seen in figure l(b) . This fairing was completely covered by the center 
section of tank B. 
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Figvire 10 indicates the manner in which additional weights were 
mounted in the tanks to simulate a given fuel load. As previously noted^ 
no effort was made to take into accoimt the effective changes in the tip- 
tank inertial characteristics due to fuel sloshing. The tip-tank weight 
is hereinafter referred to in terms of the weight ratio p which is 
defined as the ratio of the tip-tank weight to wing-panel weight by the 
relation 




Ww ^ 


since = 10 pounds . 

Figure 11 is included for the purpose of relating p to full-scale 
fuel loads for both size tip tanks. In view of the results of refer- 
ence l4 concerning the effective weight of sloshing fluid in a tank, the 
fuel loads simulated in the present investigation by means of solid weights 
correspond to increased actual fuel loads . Althoiagh tank B had three 
times the internal volume of tank A, the full-scale fuel capacity of both 
tanks was considered to be the same . The empty weight of tank B was based 
on a further assimiption that the larger full-scale tank housed fixed 
equipment whose weight was approximately 50 percent of the weight of its 
fioll fuel load. The weight of this fixed equipment accoimts for the large 
differences in p corresponding to the tank-empty condition in figure 11. 
The minimxmi weight condition indicated by the dotted boundary in the fig- 
ure, and applying to a solid fuel load of 25 percent full, corresponds 
to the lightest weight of the flutter-arresting device for tank A with 
fin. (See table I.) 

Fuselage .- The main structural member of the fuselage was a hollow 
steel section to which the wings were mounted through steel brackets 
welded to the root of each wing spar . This center section was made hollow 
so as to permit passage of the instrument wires leading from the model 
to the control room. The fairings forming the external contours of the 
fuselage were attached to either end of the main member, as may be seen 
in figure l(b) which shows model with middle fairings removed. The elastic 
properties of the fuselage were not intended to be representative of cur- 
rent fighter designs; the fuselage was, in fact, essentially rigid with 
respect to the wings . Bearings were housed at either end of the center 
section to permit a rolling degree of freedom of 1180° . Low stiffness 
springs were housed within the fuselage to position the model horizontally 
in the tunnel when the model was given freedom to roll. Small aerodynamic 
control tabs (fig. 1(a)) were located at the rear end of the fuselage to 
provide manual control of the model in roll. The rolling degree of free- 
dom could be completely locked out so as to obtain essentially the canti- 
lever or symmetric boundary condition. The actual measured rolling inertia 
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of the model fuselage was 2.56 in-lb-sec^ which is somewhat larger than 
the scaled value of a typical fighter fuselage (the scaled rolling inertia 
of the Lockheed F-94B is approximately 1.59 in-lb-sec2) . 

Instrmentation All the frequencies in both the natural vibration 
survey anrH the flutter investigation were measured by means of resistance 
wire strain gages located at three different stations along each wing 
spar as indicated in figure 5(a) . The bending gages were mounted to the 
inside siirfaces of the overhanging flanges, eind the torsion gages were 
movinted to the vertical web members. The signals from these strain gages 
were fed into a recording oscillograph. The position of the model in 
roll could be determined by means of a slide wire-position indicator 
located at the forward end of the main center-section member of the 
fxiselage . 

In order to determine the vertical acceleration of the fxxselage dioring 
flutter, an accelerometer was located on the main center -section member 
at the juncture of the fuselage center line and the wing elastic axis . 

High-speed motion pictures were taken of the model during flutter 
from a position outside the tunnel and in line with the span of the wing. 
The only other item of model instrumentation was associated with the 
operation of the flutter -arresting device, discussed in the following 
section . 

Flutter-arresting gear .- A desirable and important feature of wind- 
tunnel flutter testing is a reliable means of preventing a destructive 
buildup of amplitude once a flutter condition is attained. The more 
common. methods used at low speeds such as restraining wires or mechanical 
means within the flow are not too satisfactory at transonic speeds . 

Several alternative methods have been suggested. The Wright Air 
Development Center (WADC) , for instance, has successfully restrained 
rudder flutter in flight by quickly changing the mass balance of the con- 
trol surface . One of the objectives in this series of tests was to eval- 
uate a flutter-arresting device involving a rapid shift of the tip -tank 
center of gravity and its possible application to flight flutter-testing 
techniques . The shift of the tip-tank center of gravity was accomplished 
by moving a piston, within a cylinder running the length of each tank, 
from its extreme rearward position to its extreme forward position or 
conversely. The piston could not be stopped at an intermediate position. 
The position of the piston was controlled by an observer or by an elec- 
tronic device that monitored the electrical outputs of one of the strain 
gages on the wings . When the oscillating stress at that station exceeded 
a predetermined level, the piston was automatically fired to its other 
position. 

Inclijded in figure 5(a) is a simplified sketch of the flutter- 
arresting gear. Forward travel of the piston was obtained by supplying 
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a.ir und.6r prsssurs through ths rosxwaxd. ©nd. of th© cylind.©r, wh©r©as 
r©aru'axd mov©m©nt was accomplish©d by applying a vacuum. El©ctrical con- 
tacts w©r© placed at each end of the cylinder in order that the extreme 
forward and rearward locations of the pistons could be recorded by an 
oscillograph. 

The mass of the cylinder, being only 0*55 pound or 6 percent of the 
tank A fiO .1 weight, imposed no weight penalties except in cases where 
low fuel loads were to be simulated. The pistons vised weighed 0 . 40 , 

0.78, and 1.07 pounds, corresponding respectively to 4 . 8 , 9-2, and 12.7 per 
cent of the tank A full weight. 


MODEL TESTS 
Vibration Survey 


Before the flutter tests were started, natural frequencies of vibra- 
tion of the model were measured with the model movinted in a horizontal 
position on a dummy sting . The flexibility of this sting was made approxi- 
mately equivalent to that of the sting in the Langley l6-foot transonic 
tunnel so that the model frequencies obtained in this survey would be 
comparable with those obtained on the sting in the tunnel.. The model was 
0xd-ted both manually and with an electrodynamic shaker attached at 
various points along the wing . 


Wind-Tvinnel Test Procedure 

Figure l(a) shows the model mounted on the sting in the Langley 
16 -foot transonic tunnel. The flutter tests were conducted at Mach num- 
bers up to 0.94 and at Reynolds nvimbers per foot up to 4.2 x 10 °. 

In conducting the flutter tests, the tunnel speed was increased 
gradually and maintained at various speeds while the model response to 
random tunnel disturbances was observed and recorded with the flutter- 
arresting piston in the forward and rearward positions of the tip tank, 
gince only two positions of the piston were possible (extreme forward 
or extreme rearward) during the approach to flutter , the piston was left 
in the position that appeared to have the least damping as observed by 
the aforementioned response of the model to the random tunnel disturbances . 

The wind-tunnel tests were terminated either when flutter was encovin- 
tered or when the test could not proceed further, because of tunnel power 
limitations or lack of effectiveness of the flutter-arresting device . 
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When flutter was ehcoxmtered, the tunnel conditions were recorded 
simultaneously with the operation of the recording oscillograph and movie 
camera. At the same time, the piston was fired to the position which 
appeared to have more damping, and it was hoped that this would result • 
in a condition that was not in a flutter region (a condition which was 
usiially obtained) . The oscillograph records and movies were generally 
taken continuously before, during, and immediately after the firing of 
the piston . 

The model was excited manually by plucking the wing tips in bending 
and torsion prior to and often after each flutter run, and records were 
made of the bending and torsional vibration frequencies . These frequen- 
cies were measured after the flutter runs to provide checks on the struc- 
tural integrity of the model. 


RESULTS AND DISCISSION 


For convenience most of the experimental and calculated results of 
this investigation are indexed in table II. Figures 12 to l8 listed in 
this table show flutter Mach nimiber plotted against e^ for both experi- 
mental and calc\ilated results. Figures 19 to 23 listed in the table show 
frequency plotted against e^. 


Experimental Results 

Table III compares some of the results of the vibration survey with 
the results (available in ref. 10) of vibration tests on a current unswept- 
wing fighter-type airplane with tip tanks . The comparisons are made 
between model frequencies obtained d\jring the vibration survey and frequen- 
cies scaled from reference 10 for the closest wing and tip-tank inertial 
conditions available for both model and airplane . The frequencies and 
tip-tank inertial conditions taken from reference 10 were scaled according 
to the scale factors listed earlier in the section entitled "Description 
of Model." The agreement shown by the limited comparisons in table III 
for first bending and first torsion indicate how well the flutter model 
may have represented a current fighter -type airplane, namely, the 
Lockheed F-94 b. The fact that the two scaled second bending frequencies 
listed for the aijrplane were considerably higher than the model second 
bending frequencies coxild be partly due to reduced fuel rolling inertia 
effects (see ref. 15 ) which may have been present in the airplane but 
were not realized in the model. 

Table IV presents the results of the flutter tests in which the tip- 
tank inertial properties were systematically vaxied so as to obtain flutter 
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for ratios of xancoupled wing bending frequency to torsion fre- 
quency neeir unity. In this table the tip-tank weight ratio p 

and moment of inertia for each flutter run are represented by a simple 
designation. In this designation the letter "a" or "b" identifies the 
tip tank used, the set of numbers following the first dash refers to p, 
and the set of numbers following the second dash gives the uncoupled fre- 
quency ratio fjij /^ai* As previously noted, the value of p was deter- 
mined from the ratio of tip-tank weight to wing-panel weight; the fre- 
quency ratio was determined from figures 2h eind 25. The phase angle indi- 
cated in table IV is the angle by which the bending strain led the torsion 
strain. 


Effect of tip -tank center of gravity .- The effect of e^ on flutter 
speed may be seen in figiures 12 to I8 for various values of p and fre- 
quency ratio The various frequency ratios correspond to differ- 

ent moments of inertia for a given value of p. Speeds where flutter was 
encountered are represented by the solid points and the open points indi- 
cate speeds reached in the tests without encountering flutter . 

With the exception of figure 16, the results presented in figures 12 
to 18 were obtained with the fuselage locked to prevent roll; thus, these 
results are for symmetric flutter, that is, flutter with both wing panels 
oscillating in phase. These figvtres, except for figure 12, apply to the 
smaller tip tank (tank A) . In general, as may be observed, flutter speeds 
are Increased for center -of -gravity positions forward of the elastic axis 
eind tend to decrease as the ratio of uncoupled wing bending frequency to 
torsion frequency approaches unity. (See, for example, fig. l4.) 

The flutter frequencies corresponding to the flutter speeds (shown 
in figs . 12 to 16) are presented in fig\rres 19 to 25 • In addition, fig- 
ures 19 to 25 show the first and second measured wing natural frequencies 
with the model on the sting, the first and second calculated coupled fre- 
quencies, and the first bending and first torsional uncoupled frequencies. 
Each figure applies to a given tip-tank weight ratio p and a given cal- 
culated frequency ratio The fact that the experimental flutter 

frequency consistently fell between the measured first bending and first 
torsional wing frequencies of the model mounted on the sting indicates 
that the wing with tip tanks fluttered mainly in a combination of the 
first bending and first torsional structural modes . This type of motion 
was observed on the model during flutter . To show this , a typical flutter 
cycle, figure 26, is reproduced from the motion pictures taiken diiring 
flutter . 
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Effect of tip-tank weight .- A compaxison of the data of figure 15 
= 0.95; = 0.59 j with the data of figure I5 = 0.92; 

|j, = 1.005^ indicates that an increase of ^ from 0«59 to 1.005 is accom- 
panied by^only a moderate decrease in flutter speeds for any given value 
of e^. However, the flutter data given in table IV are insufficient to 
show the effect at low values of |i. 

An appraisal of the flutter-arresting gear .- In figures 12 to I5 
the effect of the tip-tank center of gravity on the symmetric flutter 
boundary indicates that, in general, the forward position of the piston 
should have resulted in a flutter -free condition. The fact that this 
proved to be the case is attested to by the motion pictures and the 
oscillograph records which were taken continuously when flutter occurred 
to provide complete time histories of the flutter response and subsequent 
damping action caused by the change in position of the piston. Figure 27 
is a typical flutter oscillogram which shows clearly how the forward 
position of the piston damped the oscillations . For all conditions where 
the tip-tank center of gravity was shifted forward of the elastic axis, 
the flutter -arresting gear was effective in damping flutter. 

Despite the effectiveness of this device, the shapes of the flutter 
boundaries were such as to impose certain restrictions on the practical 
application of the variable unbalance featiire of the device . These restric- 
tions are evident in a comparison of figures I5 and I6 which show increasing 
flutter speeds for tip-tank centers of gravity rearward of certain loca- 
tions . Figure I5 shows that it was possible to traverse a small portion 
of the flutter region into a safe region without dangerous amplitudes 
building up. However, figiire 16 shows a flutter condition where this was 
not possible. With the center of gravity approximately 5-6 inches rear- 
ward of the elastic axis, firing the piston forward established a condi- 
tion also very close to the flutter boundary as indicated by the open 
point at approximately 1.4 inches rearward of the elastic axis. Since 
this shift did not change the- center of gravity to a position in the safe 
region, this test could not be safely conducted beyond a Mach number 
of 0.57* 

A fvirther restriction in the application of the arresting device to 
flight-flutter use may be seen by examining the weight and size of a full- 
scale flutter-arresting device. If a full-scale device of this type were 
reproduced according to the scale factors applicable to the flutter model, 
the flutter arresting tube would be 6 times as long, or approximately 
12 feet . The full-scale piston would be 2l6 times as heavy, or would 
weigh approximately 250 pounds for the heaviest piston. 

The cumbersome features of a mass this large moving such a distance 
in a fvill-scale tank may be offset by the time of travel of the piston 
from one end of the tube to the other. In order for the piston travel 
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to be accomplished in the same number of flutter cycles for both model 
and airplane, the full-scale piston would have to traverse the arresting 

tube length in about 2^ seconds . 

2 

It is recognized that flutter conditions can occur which may be 
insensitive to a change in mass unbalance and that \mder such conditions 
more suitable flutter-arresting devices may be those which suddenly cause 
a large increase in the damping of the structure. However, for the sym- 
metric flutter encoimtered in the present study, the shifting mass- 
unbalance principle proved very effective within the limitations discussed 
in the preceding paragraphs . 

Effect of tip-tank size .- The effect of the. size of the tip tank on 
flutter can be seen in figiore 12 by a comparison of the ciurve through the 
large symbols (pertaining to the large tank) with the curve through the 
small symbols (pertaining to the small tank), both tanks having nearly 
the same value of p. and moment of inertia (resulting in nearly the same 
frequency ratio j • The reduction in flutter speed for the large 

tip tank, which was as much as 15 percent in the region covered by the 
tests, is probably due mainly to the increased aerodynamic forces resulting 
from the larger tank, although it shoxild be noted that approximately 

li percent of the area of the half -wing was covered by the large tip tank. 

2 

(See fig. 3(a)) ‘ 

Effect of tip-tank fins . - Most of the flutter r\ms were conducted 
with a horizontal fin attached near the trailing edge of the tip tank on 
the outboard side . Some tests were conducted with the fin removed, and 
the effect of this alteration on the flutter speed may be seen by com- 
paring certain tests in table IV for both large and small tip tanks . In 
some of these cases the fin appeared to be beneficial and in other cases 
it was not. For tip tank A, for example, comparisons between tests 
and 5 for a low tip-tank weight and between tests 25 and 26 for a high 
tip-tank weight show that the use of the fin significantly increased the 
flutter speed. In contrast, however, comparisons between tests 6 and 7 
for tip tank A and between tests 40 and 4l for tank B indicate that the 
addition of fins had little effect on the flutter speed and actually 
reduced the flutter speed slightly for the large tank. 

Observations on fuselage freedoms .- The effect of freeing the fuse- 
leige to roll about its longitudinal axis is shown in figure l6. These 
data pertain to the large tip teink for p = 1.02 and an uncoupled sym- 
metric frequency ratio of l.O^J-. The solid curve through the circles cor- 
responds to the fiiselage-locked condition, whereas the squares apply to 
tests in which the rolling freedom was permitted. This additional free- 
dom produced little or no effect on the flutter speed in the region where 
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symmetric flutter was encountered. Here, as in previous figures, the 
open points represent speeds reached without encountering flutter. 

Previous experiences with antisymmetric flutter (see, for example, ref. 15) 
have indicated that this type of flutter can be critical for forward tip- 
tank centers of gravity. In that case, firing the piston rearward would 
be necessary to move away from the flutter boundary and arrest flutter. 

For these tests, however, antisymmetric flutter was not encountered over 
the range of tank weights and frequency ratios covered. The failure to 
locate this boundary may be attributed, first, to the fact that the 
uncoupled first antisymmetric bending to first torsional frequency ratio 
was 5-56, whereas the first symmetric frequency ratio was 1.04; and, 
second, to the inability of the control surfaces on the tail and the soft 
torsional spring inside the fuselage to keep the model positioned in a 
near-horizontal attitude sufficiently well to permit satisfactory tests 
at Mach numbers higher than O.80 for the fuselage free-to-roll condition. 

The term symmetric flutter as used in this paper applies to flutter 
where the wings were moving in phase and was obtained when the fuselage 
was locked to prevent roll. However, from observation and motion pictures 
made during the tests, together with oscillograms of a strain gage on 
the sting and an accelerometer on the model center line at the wing elastic 
axis, it is clear that vertical motion of the fuselage was present. This 
motion was particularly noticeable for more violent flutter responses in 
which the entire model and the sting appeared to be pitching about some 
axis in the sting support system. Subsequent static and dynamic studies 
of the sting support system revealed that the sting has considerable 
flexibility and that an axis of rotation exists in part of the supporting 
structure just beneath the tionnel floor. Some Information on the elastic 
and inertial properties of the sting support system is given in the appen- 
dix. Furthennore, the fimdamental natural frequency of the entire sting 
model system appears to be well within the range of the model frequencies . 
These dynamic characteristics of the sting mount should be kept in mind 
when interpreting the symmetric flutter results obtained in this 
investigation . 

Remarks on model flutter results and related full-scale flight 
flutter experience .- Iftider certain special tip-tank inertial loading con- 
ditions involving the use of lead weights in the tip tanks, flutter was 
experienced on the Lockheed F-8OA airplane in flight tests conducted 
by WADC. Some uncertainties existed in the actual inertial loading con- 
ditions, particularly with regard to the distribution of fuel within the 
wing. Exploring such conditions on the model would have necessitated 
removal of the flutter-arresting device and was, therefore, postponed 
until the end of the experimental program. However, an attempt to repro- 
duce the flight flutter point with the model would have had uncertain 
value since the inertial properties of the model wing did not duplicate 
the scaled values of the airplane wing (figs. 4 to 6) . (it is of interest 
to note in passing, that, at least for vibration and flutter studies, the 
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wing of the Lockheed F-80A may be considered identical to that of the 
Lockheed F-9^B.) Eiirthermore , the effect, of the sting flexibility as an 
additional degree of freedom woiold make comparisons between model and 
full-scale results difficult. For these reasons, the conducting of the 
postponed test without the flutter -arresting device was not considered 
advisable . Despite the important differences that existed between the 
model and flight flutter conditions, some results have been chosen for 
conditions which, although far removed from those of the airplane, are 
as near as could be approached in the experimental program to those that 
are believed to have existed at the time the full-scale airplane fluttered. 
These results are presented in figvire 28 in terms of Mach number at 
flutter as a fimction of the ratio of uncoupled wing bending frequency to 
torsion frequency and show that the model fluttered in the same Mach mmi- 
ber range for comparable frequency ratios . Better agreement might be 
expected if, as indicated previovtsly, more accxirate inertial properties 
of the airplane had been reproduced in the model and if the dynamic pro- 
perties of the sting support system could have been adjusted to represent 
the body degrees of freedom existing in the airplane. 


Calculations and Correlation of Experiment and Theory 

Methods of calculations . - The calculated flutter speeds and flutter 
frequencies presented in this paper were determined by application of a 
Rayleigh-Ritz type of analysis in which a combination of uncoupled vibra- 
tion modes was employed to represent the flutter mode . By use of the 
measured wing inertial and stiffness properties given in figures 4 to 8, 
these modes were computed by means of the iterative method presented in 
reference 18. The uncoupled cantilever first bending and first torsional 
frequencies computed by this method are shown as functions of p, and ^ot 

in figures 24 and 25, respectively. The corresponding mode shapes are 
shown in cross -plot form in figures 29 and 50. The uncoupled cantilever 
second bending and second torsional mode shapes and frequencies are given 
in figure 51* Figure 52 shows the xmcoupled antisymmetric first bending . 
mode shapes and frequencies . 

In addition to approximating the flutter mode by means of a limited 
nvimber of uncoupled modes, the flutter calculations also involved the 
following simplifying assijmptions ; 

(1) The oscillating aerodynamic forces and moments acting on the 
wing were those derived by Theodorsen in reference 16 on the basis of 
two-dimensional incompressible potential-flow theory (except for one 
case involving compressibility noted later) . The effect of wing taper 
was accounted for in two different ways, one of which amounted to 
. weighting or grading the aerodynamic forces and moments according to the 
taper ratio as recommended in reference I9. The station at three -fourths 
of the span was arbitrarily chosen as a reference station and the 
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F + iG peirts of the aerodynamic coefficients for k at this reference 
station were held constant' along the span, whereas the remaining parts 
of the aerodynamic coefficients were allowed to vary with span according 
to the semichord ratio b/br* This method is for convenience identified 
herein by the term "graded coefficients." In the other method, labeled 
"constant coefficients," the reference station was also at three-fourths 
of the span but all parts of the aerodynamic coefficients were held inde- 
pendent of span for k at this station. 

(2) Aerodynamic forces and moments on the tip tank were not included. 

(5) The effect of structural damping in bending or torsion was taken 
into consideration in the manner recommended in reference 20; that is, 
it was assumed that g^^^ = = Constant. For most of the calculations 

this constant was taken as zero; but for a number of cases, a value of 
0.035 was used. This value was chosen on the basis of the experimental 
vibration' data and has already been mentioned in the section on scale 
factors . 

With these assumptions included in the analysis, flutter calculations 
were performed for three different wing -root boundary conditions, two of 
which were symmetric and the third antisymmetric to correspond, respec- 
tively, to the locked euad unlocked fuselage in the flutter experiments . 

The first of these boundary conditions was the ideal cantilever -wing 
condition, \diereas the second condition was that the pitching motion of 
the sting was approximated by means of a simple spring-inertia system 
whose properties are given in the appendix and was introduced into the 
flutter calculations as an additional elastic degree of freedom. For 
the antisymmetric condition, corresponding to the unlocked fuselage, the 
calculations included the first antisymmetric wing bending modes given 
in figure 52 for two different fuselage and tip-tank rolling moments of 
Inertia. The entire program of flutter calculations has been classified, 
according to the degrees of freedom used, by the code designation given 
in the following table in order to facilitate Identification of important 
features in specific applications of the analysis : 
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Identification 

Boundary 

Aerodynamic assumptions 

Degrees of freedom 

of analysis 

condition 

for wing taper 

Fuselage 

Wing 

A-1 


Graded coefficients 


First bending, first torsion 

A-2 


Constant coefficients 


First bending, first torsion 

A-3 

Symmetric 

(cantilever) 

Constant coefficients 

None 

J First bending, first torsion 
[Second bending, second torsion 

A-k 


Constant coefficients 
( con^ressible ) 


First bending, first torsion 

B-1 

Symmetric 
(sting flexible) 

Constant coefficients 

Fuselage 

pitch 

First bending, first torsion 

C-1 


Constant coefficients 


1 First asymmetric bending, 
[First torsion 




Fuselage 

L 

C-2 

i 

Antisymmetric 

Constant coefficients 

free to roll 

First asymmetric bending 
including tip-tank 
rolling inertia. 

First torsion 


For convenience in comparing the calculated ajid experimental results, 
the analytical identifications given in this table are repeated in 
table II . In studying the effects of structural damping as stipulated 
by assimiption (5) given previously, only analysis A-1 was used. 

The calculated flutter speeds and flutter frequencies are given in 
table V, which is divided into six parts, each part being identified 
according to the code given in the preceding paragraph. Tables V(a), 

V(b), and V(c) apply to the ideal cantilever-wing boxmdary condition. 
Calculations using approximations to the flexibility in the model movinting 
system are given in table V(d), and antisymmetric flutter solutions are 
given in tables V(e) and V(f) . 

Re stilts of analyses A-1 and A-2 .- Comparison of the results given 
in table V(a) for analysis A-1 (for zero damping) with those given in 
table V(b) for analysis A-2 for the same tip-tank inertial conditions 
indicates that the use of constant aerodynamic coefficients yielded flut- 
ter speeds that were as much as 30 percent higher than those obtained 
using graded coefficients. This trend toward higher flutter speeds due 
to the former and simpler application of two-dimensional aerodynamic 
coefficients was also found by Jahn and Buxton in reference 21, in which 
flutter speeds calculated by using various constant-coefficient approxi- 
mations were appreciably higher than those calculated by using a graded- 
coefficient representation that essentially involved integration of both 
F + iG and k along the span. The calculations in this reference were 
performed for a conventional bare wing of considerably lower frequency 
ratio and somewhat higher taper ratio than the frequency and taper ratios 
of the present conf igviration . 
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The results of analyses A-1 and A-2 without damping are compared, on 
the basis of Mach nvunber, with experimental flutter results in figures 15 
to l8, except for figure 17 in which the effect of damping is also shown. 
Although for some cases the comparison suffers for lack of a more complete 
fix on the experimental flutter boundaries, notably in figures 17 and l8, 
it is nonetheless evident that the disagreement between experiment and 
calculation was enlarged as the ratio of uncoupled bending frequency to 
torsion frequency approached unity. Furthermore, there were consistent 
tendencies for the calculations to be excessively conservative when the 
tip-tank centers of gravity were within 1 inch to the rear of the elastic 
axis, and for experimental and calculated flutter speeds to converge on 
one another for tip-tank centers of gravity further to the rear of the 
elastic axis . As was pointed out in reference 11, this behavior may be 
due to the failure to include some estimated tip-tank aerodynamic forces 
and moments in the calculations. (See assumption (2) in the preceding 
section.) Figiires 12 and l6 show that flutter speeds calculated by using 
constant coefficients were in better agreement with experiment than were 
those calculated by using graded coefficients . 

Effects of structural damping .- Table V(a) and figure 17 show some 
effect of introducing a damping coefficient of 0.055 into the simple two- 
mode flutter calculations involving graded aerodynamic coefficients 
(analysis A-l) . The increases in flutter speeds, due to damping, were 
generally greater for tip-tank centers of gravity within 1 inch to the 
rear of the elastic axis than for tip-tank centers of gravity further 
to the rear* Case 1*0 for e^ = 0 shows that the flutter-speed increase 

was very abrupt . Moreover, it can be noted in this table that the 
increase in flutter speed also tended to be greater as the uncoupled 
frequency ratio approached and exceeded \mity. (Compare, for 

example, case 54 with case 47 and case 55 with case 49 for approximately 
the same values of n and e^.) This trend is in general agreement with 

that found in the extensive theoretical flutter studies of reference 22. 

The effects of this small damping coefficient materially improved the 
agreement between the calculated and experimental results, particularly 
for cases where the calculations without damping yielded excessively 
conservative answers . 

Effects of higher \mcoupled modes and compressible aerodynamic coef- 
ficients.- Table V(c) corresponding to analysis A-5 shows that the intro- 
duction of various combinations of higher uncoupled modes into analysis A-2 
produced a negligible effect on flutter speed and flutter frequency for 
the case considered, namely case 64 . On the basis of the results of ref- 
erence 5, it is reasonable to expect a greater influence of higher modes 
for conditions of greater mass unbalance (corresponding to 'high e-^ val- 
ues in the present investigation) than the mass \inbalance of case 64 . 
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Although, as previously noted, most of the flutter calciilations 
presented in this paper were based on incompressible -flow theory for the 
aerodynamic approximations, the effects of compressibility were explored 
by means of analysis A-4, which was based on constant coefficients cor- 
responding to a Mach number of 0 .7 as given in reference 25 • The ,cal- 
cvilations made by vising this analysis were also applied to case 6k, ».nf^ 
the results, listed for case 65 of table V(b) , indicate that the flutter 
speed was not appreciably affected here by compressibility. As an approxi- 
mate check on this effect, examination of figure 1 in reference 25 showed 
that the k and frequency ratio of the model for this calculation fell 
approximately in a region where the same flutter speeds could be obtained 
with either compressible (M = O.7) or incompressible (M = 0) coefficients. 

The case chosen for these studies was one for which there existed 
a wide disagreement between theory and experiment and it is apparent that 
no Improvement was obtained by Introducing either higher uncoupled modes 
or compressible aerodynamic coefficients . 

Note on coupled modes of vibration .- Along with the other frequencies 
in figures 19, 20, and 25 are shown also the first two coupled cantilever 
frequencies which were calculated from the flutter stability determinants 
with the aerodynamic terms omitted. These computations were performed 
for the purpose of estimating how much of the discrepancy between measured 
frequencies and calculated uncoupled frequencies was due to coupling in 
the system. As may be observed, coupling constitutes a relatively small 
part of this discrepancy, especially for the first mode. The remaining 
gap, of as much as 50 percent, between measured and calculated coupled 
frequencies is probably due in large measure to the fact that the measured 
frequencies include a strong influence of the sting support system. 

Effect of sting flexibility .- The effect of flexible pitching in 

the sting mount, as simulated by means of analysis B-1, is shown in 

table V(d) and figure 53- Case 6k was used in this study which included 

a range of sting frequencies fg and equivalent sting-fuselage pitching 

moments of Inertia If g • The calculated flutter speed is seen to be 

sensitive to frequency fg, but relatively insensitive to inertia If „ • 

s 


In order to compare the results of analysis B-1 with experiment, 
estimates of the parameters (ifg = 26.4 ft-lb-sec^ and fs = 16.8 cps) 

are made in the appendix euid are shown as dashed lines on figure 53* By 
using the estimated values of the parameters, it can be observed that 
there is very little effect of sting flexibility on flutter speed for 
the case considered. However, the large effect of frequency fg indi- 
cates that sting flexibility may be important when the flutter frequency 
is close to the sting frequency. 
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Antisymmetric boundary condition .- The results of the antisymmetric 
flutter calculations are given in tables V(e) and V(f) for analyses C-1 
and C-2, respectively. These calc\ilations Involved the first uncoupled 
antisymmetric bending-mode shapes and frequencies presented in figure 52 . 
The fact that the flutter speeds and flutter frequencies are consldferably 
higher than those obtained from cantilever analysis A-2 for the same tip- 
tank weight and moment of inertia is due to the greatly increased frequency 
ratios based on the antisymmetric bending frequency . An overall pictiore 
of the effect of frequency ratio in relation to tip-tank center -of -gravity 
location is shown for three frequency ratios in figure 54 and corresponds 
to analysis C-1. As may be seen, the frequency-ratio effect was consid- 
erably more pronounced for forward tip-tank centers of gravity than for 
rearward tip-tank centers of gravity. Notice also the reversed effect 
of this parameter on flutter speed as the tip-tank center of gravity is 
moved from a position forward of the elastic axis to a position rearward 
of the elastic axis . The value of Ip used in analysis C-1 was estimated 

during the design of the model. The higher value of If^ used in analy- 
sis C-2 was based on measured data obtained after completion of the model. 
However, the differences in flutter speeds between tables V(e) and V(f) 
are attributed more to the effect of I-tj. than of Ifj.* This contention 

is based on the existence of a small difference between the antisymmetric 
bending frequencies due to the two different values of If^., as compared 

with the larger difference caused by modifying the tip bending-moment 
boundary condition from zero to a finite value determined from the esti- 
mated tip-tank rolling Inertia about the wing tip. (See fig. 52.) How- 
ever, as comparison between tables V(e) and V(f) shows, the effect of 
this parameter on flutter speed was small. 

The theoretical antisymmetric flutter trend plotted in figure l6 
(for analysis C-2) serves to show the location of the experimental sym- 
metric flutter region. As is implied in figin:e 54, this region might well 
have fallen within the speed range covered by the experiments, particularly 
for forward tip-tank centers of gravity, had the ratio of lancoupled anti- 
symmetric bending frequency to torsion frequency been closer to unity. 
Preliminary analytical work during the construction of the model indi- 
cated that this condition might be realized for tip-tank weights 
approaching an empty condition, a condition which was not attainable in 
this investigation . 


CONCLUSIONS 


An extensive investigation has been made of flutter of a true-speed 
dynamically scaled model representative of current unswept-wing fighter-, 
type airplanes with tip tanks . The model was one-sixth the size of a 
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hypothetical fighter prototype having an aspect ratio of 6 and was of 
spar-balsa segment-type construction with a flutter-arresting device in 
each tip tank. The following conclusions can be drawn from the results 
of this investigation: 

1. A true-speed dynamically scaled model can be built to satisfy 

a given set of elastic and inertial specifications for flutter tests at 
high subsonic speeds . 

2. In general, symmetric flutter speeds were increased for tip-tank 
center -of -gravity positions forward of the elastic axis and tended to 
decrease as the ratio of uncoupled wing bending frequency to torsion 
frequency Ihi^^’a! approached unity. 

5. Increasing the volume of the tip tanks by a factor of 3 tended 
to decrease the flutter speed. The installation of a horizontal fin 
at the rearward end of the tank may have had a beneficial effect on the 
flutter speed for the small tank, but may have had a slightly detrimental 
effect for the larger tanks . 

4 . Symmetric flexure -tors ion flutter was effectively arrested when 
the tip-tank center of gravity was quickly shifted forward of the elastic 
axis . For cases where the center of gravity could not be shifted this 
far forward, the effectiveness was not as positive and was in some cases 
negative. In applications to flight flutter use, the weight, size, and 
time factors involved in the effectiveness of this device must be 
considered . 

5 . Calculated flutter speeds based on assumptions involving struc- 
tural damping, higher structirral modes, compressible aerodynamic coeffi- 
cients, and sting -support flexibility (and not inclviding aerodynamic 
forces and moments on the tip tanks) were, in general, excessively con- 
servative for tip-tank centers of gravity near the elastic axis . How- 
ever, the introduction of structural damping improved the agreement 
between calculation and experiment and resulted in greater increases in 
calculated flutter speeds, especially for ratios of uncoupled wing bending 
frequency to torsion frequency near unity. The effects of sting-support 
flexibility showed that this modification to the calculation may be impor- 
tant where the natural frequency of the sting system was near the flutter 
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frequency, and further work on these effects is in order. In addition, 
this sting flexibility should be kept in mind in interpreting the experi- 
mental flutter results . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 29, 195^* 
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APPENDIX 

CHARACTERISTICS OF MODEL MOUNTING SISTM 
Description 


The model was mounted in the Langley l6-foot transonic timnel in 
the manner shown in figure l(a) . This mounting system essentially con- 
sisted of two main parts, the. sting and the sting-support structure. 

The junctijre between these two parts was 95 inches rearward of the elastic 
axis along the tunnel center line, and the actual connection was made by 
mpgi ns of bolts approximately 1 inch in diameter located ^ound the periph- 
ery. The sting had a joint whose center line was located 60 inches rear- 
wsird of the elastic axis along the tunnel center line. At the free or 
upstream end of the sting a sting extension was mounted which was part 
of the main structural member of the model fuselage. 

The sting support structvire was made up of a massive streamlined 
strut mounted to a circular -shaped track that was part of the mechanism 
used to provide a change in angle of attack on models in the tvinnel. The 
top of the strut to which the base of the sting was connected was cylin- 
drical in shape with the longitudinal axis of the cylinder lying along 
the tunnel center line . 


Stiffness Distribution 

Figure 35 gives the stiffness distribution along the length of the 
sting . The variation shown is based on both calculation and experiment , 
the experimental, values being obtained from measurements • of the bending 
slope with the sting subjected to vertical loads applied to the main 
structural member of the f\iselage at the elastic axis . Along the solid 
portion of the curve the calculated and experimental values were- in close 
agreement. The dotted portions of the curve apply to the various joints 
and connections in the system and represent the most reasonable stiffness 
values • corresponding to the actiial measured slopes over these discontinui- 
ties. The experiments further revealed that the sting support structure 
was, for all practical purposes, rigid in the vertical or pitch direction, 
but that there was a finite stiffness contribution from the strut mount 
on the circular track. An estimation of this stiffness in terms of a 

spring constant for the pitching degree of freedom is 596 x 10^ in-lb/radian 
based on the approximate location of the actioal axis of rotation in the 
base of the sting support structure . 
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Weight Distribution 

Figvrre 56 shows the variation of weight per unit length with the 
distance along the tunnel center line. This distribution was calculated 
from the specifications of the sting. 


Fuselage Inertial Properties 

The weight and pitching moment of inertia of the fuselage, including 
the main structural center section, wing spar root brackets (see fig. 5(a)), 
nose, center, and rearward fairings with control-surface assembly, are 
given as follows : 


Fuselage weight, lb I 5 Y 

Mass pitching moment of Inertia of fuselage about 

elastic axis, In-lb-sec^ 75.3 

Fuselage center of gravity, rearward of elastic 

axis, in 5.92 


Simplified Representation of Sting Flexibility 

A scheme is advanced herewith for approximating the actual elastic 
behavior of the sting mount by means of an equivalent one -degree -of - 
freedom system oscillating in pitch about a point S in the manner shown 
in the following sketch: 
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The equivalent elastic properties consisted of a spring constant kg 
and an effective rotation arm rs . These parameters were determined 
from static loadings at the elastic axis together with deflection-curve 
measurements along the sting. As indicated in the foregoing sketch, the 
slope of the deflection curve at the elastic axis established the point S 
as being the effective axis of rotation and Tg as the effective rota- 
tion arm. From these measurements the following values were found for 
these two parameters : 


kg = lt015 Ib/in. 

and 


rg = k2 in. 


In addition. 


If. = 26. k ft-lb-sec2 
Is 

which is based on the measiired fuselage inertial properties given pre- 
viously and includes the effective sting mass corresponding to 
Tg = k2 incihes • 

On the basis of the foregoing structural and inertial properties^ 
the estimated frequency of the fuselage -equivalent sting combination is 


^s 


16.8 cps 

2. ^24lf^ 
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TABLE I.- GEOMETRICAL AND INERTIAL PROPERTIES OF MODEL TIP TANKS 


(a) Small tip tank: tank A (scaled 230-gallon tip tank) 


[Length, 30 -5 in-; maximum diameter, 4.24 in.; 
span of fin (from tank center line), 4-38 in.; 
exposed area of fin, 11.2 in.^ 


Piston 

weight, 

lb 


Without 

fin 


With 

fin 

Wt, 

lb 

in- 

in-lb -sec^ 

Wt. 

lb 

in. 

I(Xt> 

in-lb -sec^ 

0.40 

2.70 

1.46 

0.280 

2.78 

1.96 

0.352 

•78 

3.08 

2.49 

.378 

3.16 

2.92 

.430 

1.07 

5.57 

3.22 

.466 

3.45 

3.59 

.520 


(b) Large tip tank: tank B 

[Length, 43-7 in.; maximum diameter, 6.12 in.; 
span of fin (from tank center line), 6-31 in.; 
exposed area of fin, 23-9 in.^ 


Piston 

weight, 

lb 


Without 

fin 

With fin 

Wt, 

lb 

®t^ 

in. 

in-lb-sec^ 

Wt> 

lb 

et^ 

in. 

^at^ 

in-lb-sec^ 

0.40 

4.77 

1.22 

0.698 

5.00 

1.92 

0.803 

.78 ■ 

5.15 

1.75 

.680 

5.58 

2.48 

.901 

1.07 

5.44 

2.22 

.769 

5.67 

2.92 

.990 
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TABLE II.- INDEX TO BESUUPS 


Table 

Description 

III 

Experimental - comparable natural frequencies of model and a typical 
, full-scale fighter configuration 

IV 

Experimental - resiilts of model flutter program 

V 

Analytical - results of flutter calculations 


Figure 

Tank 

size 



In. 

Wing-root 

boundary 

condition 

Analysis 




Ordinate , 

M; abscissa. 

«t 


12 

A and B 

0.60 

0.92 

- 4.9 to 1.7 


None 

15 


•59 

• 84 , .93 

-4.9 to 3-5 


A- 1 , A -2 

14 

A 

.82 

•90, .96, 1.02 

-5.1 to 4.0 

Fuselage 

clamped 

A -1 

15 


1.005 

•92 

- 2.4 to 4.0 

A -1 




1.04 

- 2.4 to 4.0 


A- 1 , A -2 

16 


1.02 

3.36 

- 2.6 to 1.6 

Fuselage 

None 


B 


3-36 

0 to 2.4 

free to 
roll 

C -2 

17 


.895 

•98 

-3.8 to 4.0 

Fuselage 

clan 5 )ed 

A -1 

18 

i 

1.29 

1.01, 1.09 

- 2.1 to 4.0 

A -1 




Ordinate , 

f; abscissa. 

et 





.84 

-1.0 to 2.8 


A -2 

19 


.595 

.84 

- 2.0 to 3.5 


First two coupled 
freqiiiencies 


A 


•93 

- 4.8 to 1.7 


A -2 

20 


.59 

•93 

-1.25 to 1.6 

Fuselage 
- clamped 

First two coupled 
frequencies 

21 


.82 

1.02 

-.9 to 2.5 

A -1 

22 


1.005 

•92 

-2.3 to 4.0 


A -1 




1.04 

- 2.4 to 4.0 


A -2 

23 

B 

1.02 

1.04 

0 to 3.2 


First two coupled 
frequencies 


CONFIDENTIAL 







TABLE III . - COMPARISON OF MODEL VIBRATION DATA WITH SCALED 


54 


CONFIDENTIAL 


NACA RM L54I19 


E3 

o 

•H 

-P +> 
m -H 
(D ^ 

o 


oj m 
ft 

ft a 


H c /3 
•p ft 
ft o 


H CQ 

rQ ft 
ft c; 


cvj 

a 

<D 

•\ to 

I 

a 


p a 

QJ *H 




oJ cd 
q tio 
• I 


H < H < 


I 00 

I • 

I VO 

I m 


VO 

d 


VO 


rovoo 

OJ 

K^ to 


ir\ OJ 


O o 

CVJ OJ 


ovoo 

H rH 


00 

moj 
to to 


VO H 

^5 


VO 

ov 


q H 
to to 


I lO 

I 

I 


!8 
I H 


H pq 
H 


I CVJ 

I • 
I 

I VO 


00 

to to 

00 O 
rH 


to lO 


OV H VO 

CVJ to toto 


ir\^ VO ON 


00 

iH rH 


CVJ to 
VO ^ 


00 o 
O 

to 


-4-10 O QV ^ Q> 

t^OO CO 00 


to to to to 


I O 
1 Q 

i P 


• - • 

• • 

* 

• 

■p 

• (U 

i 

• OJ 

i 

• 

0 

§ 

p 

0 

p 

o 

H H 

H H 

p 

p 

0 ft 

0 ft 

0 

ft 

s 

P M 

P U 

p 

u 


O P 

O P 

o 

p 


s < 

s < 

s c 

H 


l 30 

P 

ft 

(U 

0 

CO 


0 

rH 

p 


ir\ 

p 

CQ 

0 

P 

o 

u 

ft 

cd 

td 

Q 


t 30 

•H 

ft 

0 

0 

CO 


p 

p 

0 

c 

o 

p 

0 

p 




0 

p 

O 

s 


H 

H 


B 0 P 
0 bD 


pq 


CONFIDENTIAL 



NACA RM Ta^kllS 


CONFIDEWTIAL 


35 



CONFIDENTIAL 


2.06 



TABLE IV.- RESULTS OF EXPER3MENTAL FLUTTER PBOGRAM - Continued 
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violent flutter; symmetric mode Involving fuselage translation 

No flutter 

Flutter 

Flutter 

Visual evidence of flutter; piston In forward position appeared very effective in arreatlng flutter 

No flutter; torsional oscillations appeared stronger and less damped with piston forward than 
with piston rearward 

Very close to flutter; piston appeared Ineffective In arresting flutter 

No flutter within range of teat; fitting for left wing taidc loosened during test 

Flutter 

Flutter 

Flutter 

Flutter 

Flutter 

Flutter 

Flutter 

Flutter In bursts 

No flutter within range of available tunnel power; action of piston very effective in preventing 
flutter 

Flutter; action of piston very effective in preventing flutter 

Flutter; piston effective In stopping flutter 

Flutter 
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Remarks 

Flutter 

No flutter; wing shoved strong bending tendencies 

Flutter; strong torsion response 

Flutter 

Flutter; amplitude of oflclUations increased slowly prior 
to flutter 

Flutter 

No flutter; response of model observed during test 

appeared as close to flutter with piston forward as with 
piston rearward 

Test discontinued; no flutter observed 

No flutter; xxredominant bending response 

Sustained flutter; predominantly torsion; piston very 
effective 

No flutter; model free to roll 

No flutter; model free to roll; very unstable In roll; 
records indicate model less damped with piston In 
forward position near top Mach nximbera 

Flutter; model free to roll; piston forward effective in 
stopping flutter 

Flutter in bursts; model free to roll; strong torsion 
response 

Phase 
angle , 

to 
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TABLE V.- RESULTS OP PLOTTER CALCULATIONS 


(a) Cantilever -wing boundary condition, graded coefficients, 
first bending and first torsion (analysis A- 1 ) 


Case 

et^ 

in. 



0 

Shi 

= gai= 0.035 

Vf, 

fps 

ff. 

cps 

l/kr 

Vf. 

fps 

ff. 

cps 

l/kp 



A 

-.565- *91 configuration 



1 

3.925 

699 

20.8 

13.09 





A-.59-.93 configuration 

2 

-1.16 







5 

0 

501 

22.7 

8.60 




4 

.24 

445 

22.6 

7.66 

562.5 

22.25 

9.85 

5 

•50 

46 l 

22.5 

8.00 

523 

22.2 

9.18 

6 

1.00 

519 

22.2 

9.12 

547 

22.1 

9.63 

7 

1-59 

558 

21.8 

9.96 

583.5 

21.95 

10.35 

A-. 59-.84 configuration 

8 

0 

993 

24.1 

16.05 




9 

.50 

774 

23.7 

12.74 

849 

23.2 

i 4.28 

10 

1.00 

748 

23.3 

12.51 

774 

23.0 

13.1 

11 

1-75 

761 

22.6 

15.13 

785 

22.7 

15.45 

12 

2.64 

785 

21.8 

14 .04 

851 

22.3 

14.54 

15 

3.50 

803.5 

21.1 

14.84 

886 

21.9 

15.75 




L-.8i5-.9O configiiration 



l 4 

0 

655 

19-9 

12.82 




15 

1.00 

561 

19.4 

11.26 

582 

19.3 

11.74 

16 

2.00 

615 

18.6 

12.9 

658.5 

19.0 

15.51 

IT 

3.00 

652 

17-9 

14.2 

729 

18.5 

15.51 

18 

3.695 

667 

17.4 

14.95 

777 

18.5 

16.56 

19 

4.00 

672 

17.2 

15.23 

781 

18.0 

16.9 




l-.825-.96 configuration 



20 

0 

329 

18.9 

6.78 




21 

.50 

367 

18.8 

7.60 

422 

18.65 

8.80 

22 

1.02 

442 

18.5 

9.32 

475.5 

18.6 

9.98 

23 

1-97 

516 

17-9 

11.25 

574 

18.3 

12.23 

24 

3.00 

563 

17.2 

12.75 

656 

17.9 

14.27 

25 

4.00 

592 

16.6 

13.9 
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TABLE V.- RESULTS OF FLUTTER CALCULATIONS - Continued 


(a) Cantilever -wing boundary condition, graded coefficients, first bending 
and first torsion (analysis A-1) - Concluded 




Vf, 

fps 

cps 

H 

Vf, 

fps 


A-. 82 -1,02 configuration 


5-1^9 

6.85 

1*05.5 

18.0 

8.72 

8.40 

1*51* 

18.05 

9.80 

9.82 

11.11 

11.8 

12.98 





A-1 .005“ -92 configuration 


B-.895-'98 configuration 


0 

H 0 

1 

217.5 

18.0 

4.705 

.69 

376 

17.8 

8.22 

1.50 

472 

17.3 

10.64 

2.50 

536 

16.7 

12.5 

i*.oo 

585 

15.8 

14.45 


45 

..27 

46 

.0135 


.65 

48 

.94 

49 

1.96 

50 

3.58 

51 

3.99 



52 

-•35 . 

53 

.50 

54 

1.50 

55 

2.50 

56 

4.00 


57 

.99 

58 

1.50 

59 

2.50 

60 

4.00 


B-l.02-l.04 configuration 


B-1.29-l*01 configuration 




B-1.29-1‘09 configuration 



9-54 

501 

14.15 

10.33 


11.92 



13.77 
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TABLE V.- EESULTS OF FLUTTER CALCULATIONS - Continued 


(b) Cantilever-wing boundary condition, constant coefficients 
first bending and first torsion (analysis A-2) 


Case 

^t^ 

in. 

Vf, 

fps 

ff, 

cps 

l/kr 

A- -565- *91 configuration 

61 

5.925 

762 

20.9 

14.25 

A-.59-*95 configuration 

62 

-1.16 




63 

0 

664 

22.7 

11.39 

64 

.24 

486 

22.6 

8.39 

65a 

.24 

488 

22.2 

8.67 

■ 66 

•50 

491 

22.6 

8.50 

67 

1.00 

524 

22.8 

8.97 

68 

1.59 

6o4 



21.9 

10.8 

A- -59- *8^ configuration 

69 

-2.00 




70 

0 

999 

24.1 

16.16 

71 

.50 

831 

23.7 

13.66 

72 

1.00 

813 

23.3 

13.61 

73 

1.75 

819 

22.6 

l4.l4 

74 

2.64 

855 

21.8 

15.3 

75 

3.50 

8^ 

21.1 

L J 

16.45 

B-l .02-1.04 configuration 

76 

-.27 




77 

.0135 

720 

16.2 

17.35 

78 

.50 

374 

16.1 

9.05 

79 

.94 

407 

15.9 

9.95 

80 

1.50 

442 

15-9 

10.86 

81 

1.96 

475 

15.6 

11.86 

82 

2.75 

511 

15.2 

13.08 

83 

3.22 

532 

15.0 

13.82 

84 

3.99 

556 

14.7 

14.71 


^Computed using Mach nimber 0.7 coefficients 
(analysis A-4) 
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TABLE V.- RESULTS OF FLUTTER CALCULATIONS Continued 


(c) Cantilever wing boundary condition, constant coefficients, first bending, 
first torsion, second bending, and second torsion 

[(Analysis A-3) ; configuration A-.59-*95; = 0.2if in^ 


Case 

Identification of inodes used 

Vf, 

fps 

ff, 

cps 

l/kr 

85 

First bending, first torsion, 
second bending, second torsion 

479 

22.7 

8.25 

86 

First bending, first torsion, 
second bending 

476 

22.6 

8.20 

8? 

First bending, first torsion, 
second torsion 

476 

22.6 

8.20 

64 

First bending, first torsion 

486 

22.6 ' 

8.59 


(d) Symmetric (sting-fuselage pitching) boundary condition, constant 
coefficients, flexible sting pitching, first cantilever wing 
bending, and first cantilever wing torsion 

[(Analysis B-1) ; configuration A-.59-*95; ©t = 0.24 in.J 


Case 

ft-lb-sec2 

Vf, 

fps 

cps 

lAr 

fg = 14.5 cps 

88 

25 

465 

25.0 

7.87 

89 

39.8 

431 

23.0 

7.29 

90 

50 

427. 

23.3 

7.13 

91 

TO 

435 

22.75 

7.45 

92 

100 

447 

22.8 

7.63 

fg = 16.7 cps 

95 

10 

518 

22.7 

8.90 

94 

25, 

509 

23.05 

8.61 

95 

40 

463 

25 .0 

7.82 

96 

50 

428 

22.9 

7.28 

97 

TO 

447 

22.8 

7.64 

98 

100 

426 

22.9 

7.26 

fg = 20 cps 

99 

10 

46l 

22.9 

7,84 

100 

25 

658 

23.1 

11.075 

101 

39.8 

536 

23.1 

9.05 

102 

50 

528 

25.1 

8.91 

103 

TO 

469 

22.9- 

7.97 

fg = 25 cps 

104 

10 

995 

23.0 

16.85 

105 

25 

1,112 

23.3 

18.58 

106 

40 

970 

24.35 

15.52 

107 

. 50 

990 

23.1 

16.72 

108 

TO 

999 

23.0 

16.90 
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TABLE V.- RESULTS OF FLUTTER CALCULATIONS - Continued 


(e) Antisymmetric (ftiselage free to roll) boundary condition, 
constant coefficients, first asymmetric bending, and first 
cantilever wing torsion 


jAnalysis C-1; configuration B-l.02-l.04; 

If^ = 0.0605 ft-lb-sec2; I^^ = ^ 


Case 

et^ 

in. 

fps 

ff, 

cps 

1/kr 

109 

-2.50 

3,441 



110 

-2.00 

4,213 



111 

-1.00 

7,125 

51.3 


112 

-•75 

9,995 

51.0 

125.7 

113 

-.50 

2,829 

15.4 

71.7 

114 

0 

1,427 

15-1 

36.7 

115 

1.19 

1,088 

14.6 

29.2 

116 

2.575 

1,027 

14.0 

28.6 


TABLE V.- RESULTS OF FLUTTER CALCULATIONS - Concluded 


(f) Antisymmetric (fuselage free to roll) boundary condition, 
constant coefficients, first asymmetric bending, and first 
cantilever wing torsion 

jtoalysis C-2; configuration B-1. 02-1.04; If^ = 0.0984 ft-lb-sec^; 

Itj. = 0.0192 ft-lb-sec2] 


Case 

in . 

Vf, 

fps 

ff, 

cps 

1/kr 

117 

0 

1,489 

15.05 

38.6 

118 

1.19 

1,131 

l4.4 

30.6 

119 

2.575 

1,059 

13.8 

29.9 
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Figure 1.- Photographs of dynamic model. 



kk 


CONFIDENTIAL 


NACA RM L54I19 



(b) Model mounted on dumny sting for vibration sumrey. L-77^39 
Figure 1.- Concluded. 
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Figure 2.- Tip tanks used in investigation. T‘-77Vl3 
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Wing semispan, inches 

Figure 4.- Spanwise weight distribution. 
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Figxrre 5.- Variation of center of gravity and elastic axis with span. 
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Figure 7.- Variation of bending stiffness with span. 




Figure 8.- Variation of torsional stiffness with span. 
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Wing semispan, inches 

Figure 9*- Design specifications for wing spar. 


CONFIDENTIAL 


54 


CONFIDENTIAL 


NACA RM L54I19 



CONFIDENTIAL 



MCA RM L54I19 


CONFIDEHTIAL 


55 
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Fuel load, percent full 


Figvire 11.- Correlation of model tip-tank weight ratio with prototype 

fuel load. 
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Figure l4.- Variation of experimental and theoretical flutter Mach numbers 
vlth tip-tank center of gravity for three values of Tank A; 

^ = 0.82; fuselage clamped. 







Figure I5.- Variation of experimental and theoretical flutter Mach numbers 
with tip-tank center of gravity. Tank Aj ^l = I.OO5; fuselage clamped. 
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O No flutter 
• Flutter 

□ No flutter" 
■ Flutter 

Fuselage free 

1 Fuselage clamped 

■ Analysis A-2 
Analysis A-l 

■ Analysis C-2 




O 


o 

□ 


□ 


Elastic axis 


Forward 


Rearward 


Figure 16.- Variation of experimental and theoretical flutter Mach nvunbers 
with tip-teink center of gravity with fuselage clamped and fuselage free 
to roll. Teuik B; p. = 1.02} ~ 1»04. 
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Figure 17*- Variation of experimental and theoretical Mach numbers with 
tip-tank center of gravity for f^j^ /^ai - 0 .98. Tank B; n = O.895; 
fuselage clamped. 
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Figure 19.- Comparison of analytical and measured frequencies as a function 
of tip-tank center of gravity. Tank A; p = O.595; ^hi/fai = O.Slf; 

fuselage clamped. 
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Figure 20.- Con^jarison of analytical and measured frequencies as a function 
of tip-tank center of gravity. Tank A; ji = O.59; = 0.93; 

fuselage clan^jed. 
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Figure 21.- Conrpeirison of analytical and meeisured frequencies as a function 
of tip-tank center of gravity. Tank A; p = 0.82; ^hij^a.i ~ 1*02; 

fuselage clan5)ed. 
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Figure 22 .- Comparison of analytical and measured frequencies as a function 
of tip- tank center of gravity. Tank A; p = 1.005; = 0.92; 

fuselage clan^jed. 
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Figure 23-- Con5)arison of analytical and meas\ired frequencies as a function 
of tip-tank center of gravity. Tank B; p = 1.02; 

fuselage clamped. 
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Figure 2k. - Variation of first 

we 


C 


MCA RM L54I19 


CONFIDEKTIAL 


69 
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2 

, m- lb- sec 

Figure 25.- Variation of first uncoupled torsional frequency with tip- 

tank moment of inertia. 
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^ l-861;39 

Figure 26.- High-speed motion pictures showing a typical cycle of symmetric 
flutter. The direction of increasing time is from top to bottom begin- 
ning with frame 1. 
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Figure 27-- Portion of oscillogram taken at flutter showing effect of 

flutter-arresting device. 
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in-lb-sec^ 


Figure 50.- Cross plot of uncoupled cantilever 
function of tip-tank moment of 

0 


torsion-mode shape as a 
inertia. 
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”^sym 


Figure J>1.~ Second uncoupled bending and torsion-mode shapes. Tank A; 

“ 0-59j ~ ®*95* 
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(a) Effect of sting- fuselage pitching moment of inertia on flutter 

speed for four values of sting- fuselage frequency f 3 . 



(b) Effect of sting-fuselage frequency on flutter speed 

for If = 26 . k. 
s 

Figure 55 •- Correlation of experimental with theoretical flutter speeds 
including effects of sting flexibility (analysis B-l). 
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Figure 5^.- Calculated antisyimnetric flutter velocity (analysis C-l) 
against tip-tank center of gravity for various values of the 
frequency ratio. Tank B; |i = 1.02; If_ = O.O 605 . 
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Distance from wing elastic axis, inches 


Figure 35 *“ Sting stiffness distribution. 
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